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'^ ■ Abstract 

^ ■ We present an interpretation of recent cyclotron resonance experiments on 
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the two-dimensional electron gas in GaAs/AlGaAs heterostructures. We show 



f~^ ■ that the observed dependence of the resonance spectrum on Landau level 

0^ . occupancy and temperature arises from the interplay of three factors: spin- 



splitting of the cyclotron frequency; thermal population of the two spin states; 



B 

i-rt ■ and coupling of the resonances for each spin orientation by Coulomb interac- 

O ■ tions. In addition, we derive an /-sum rule which allows spin polarisation to 



be determined directly from resonance spectra. 
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Since the discovery of the quantum Hall effect, there has been intense interest in the 
properties of a two-dimensional electron gas in a strong magnetic field. Although it was the 
dramatic transport phenomena that originally drew attention to these systems, increasing 
use is being made of other experimental probes. At very low filling fraction, transport mea- 
surements are problematic and optical techniques such as photoluminescence and cyclotron 
resonance are particularly valuable. The optical spectra show many interesting features, and 
are still far from understood. 

Recently, two groups have reported intriguing cyclotron resonance measurements in GaAs 
inversion layers at low Landau level filling fraction, v 0,@]. For v > 1/6, a single resonance 
peak is observed, which shows no anomalous behaviour at any of the fractional quantum 
Hall states. Below u ^ 1/6, a second resonance peak appears, and the absorption spec- 
trum is strongly dependent on both filling fraction and temperature. These observations 
were interpreted initially as a signature of Wigner crystallization [|l| . Subsequently, it was 
demonstrated |[1|,^ that the frequency splitting of the two resonances at the lowest fill- 
ing fractions is very similar to the spin-splitting measured for cyclotron resonance in bulk 
GaAs, which arises because the effective (^-factor is energy-dependent [0]. It has since been 
suggested that the observations may indicate a spin-ordering transition, or localization of 
electrons in the extreme quantum limit 0. 

In this paper we offer a detailed theoretical interpretation of these experimental results. 
We show that the filling-factor dependence of the cyclotron resonance spectrum is not due 
to any change in the positional or magnetic correlations of the two-dimensional electron gas. 
Instead, it is the result of a cross-over from independent, spin-split resonances at low u, to a 
single mode dominated by Coulomb interactions at high z/. Moreover, we are able to account 
for the temperature dependence of the spectrum in terms of the thermal populations of the 
two spin states, calculated using the known single-particle f^-factor and without invoking 
exchange interactions. 

Kohn's theorem |Q provides a natural point of departure for theories of cyclotron line- 
shape: in a translationally invariant, one-component system, cyclotron resonance couples 



only to the centre of mass, and the resonance is unaffected by interactions or correlations. 
Spectral structure may be observed if translational invariance is broken, either by disorder 
1^ or by non-parabolicity of the conduction band . Alternatively, it may arise if there are 
carriers of more than one type, with different cyclotron frequencies, so that centre of mass 
motion does not separate from internal motion of the electron gas. This happens in strained 
Si inversion channels, in which separate pockets of the conduction band are occupied, and 
has been studied experimentally and theoretically |^ at weak magnetic field. There is a 
clear indication that the same mechanism is at work in the experiments we are concerned 
with, from the fact that the frequency splitting of the resonance at low v matches the bulk 
spin-splitting 0. 

Motivated by this, we develop in the following a theory of cyclotron resonance for a 
two-component system in the strong field limit, the two components representing spins of 
each orientation. A complete treatment of many-body correlations at finite temperature 
would clearly be very difficult. Fortunately, the experimental conditions suggest several 
simplifications. First, the filling fractions of interest are small, so exchange energies are 
negligible [0, and the electrons may be treated as distinguishable particles, each with a 
conserved spin orientation. Second, the cyclotron energy is much larger than all other 
energy scales, so that only the two spin states of the lowest Landau level are populated, 
and cyclotron resonance involves transitions only into the corresponding spin states of the 
next level. Since higher Landau levels do not contribute, the two spin orientations may each 
be treated as if they have a parabolic dispersion relation, with slightly different effective 
masses. The relative concentrations of these two types of particle are fixed simply by their 
Boltzmann weights, in terms of the Zeeman splitting of the lowest Landau level. The 
favoured spin state has the higher cyclotron frequency, so the minority spins can be thought 
of as randomly distributed 'heavy' impurities in a sea of the majority spin state. 

At low density, the natural co-ordinates are the guiding centre and orbital co-ordinates 
of the electrons 



R. = n-p^,a^ = '-^,al = 'l±Jfl, (1) 



where Pj = (pj + eAi) x z/eB, and / = Jh/eB is the magnetic length. The only non-zero 
commutators are 

[ai,al] = 6ij and [Ri,B!j]=il^6ij. (2) 

Cyclotron resonance measures the real part of the dynamical conductivity, which is propor- 
tional to iig 

1 fOO 

S{iu) = - 5:(a,(t)a](0))e'-*dt. (3) 
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In view of the low filling fraction, we expand the interactions in powers of the orbital co- 
ordinates. At leading order in the interactions, the equations of motion are 



• _ _1_^ dU_ 
eB dvj 
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where tu, is the cyclotron frequency of electron-i, U{{r^}) is the Coulomb energy, rapidly 
oscillating terms are omitted and an additive constant has been absorbed in a,. 

Thus, from Eq ^, the guiding centres drift along contours of the local potential gradient 
(independently of the orbital motion), while, from Eq |], the orbital co-ordinates perform 
coupled oscillations, with the coupling determined by the instantaneous positions of the 
guiding centres. For Coulomb interactions the coupling matrix is 

M,, = \ ^ (6) 

At the lowest experimental temperatures and densities, the system is expected to form 
a Wigner crystal, with the guiding centres located at the sites of a triangular lattice. In 
these circumstances, the coupling matrix is time independent and the resonances are normal 
modes satisfying 
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where we measure frequency, tu, relative to the cyclotron frequency, cjc, of the majority 
spins, in units of the splitting, dujc-, so that 5uJi = 0, —1 for the majority and minority spins, 
respectively. The dynamical matrix (^) is defined by the triangular lattice vectors, in units 
of the lattice constant, a = {^/S/2 n)~^/^, and / is the dimensionless interaction strength 

Sneeoa^ MuJc 

Hence we have established a model for cyclotron resonance in the Wigner crystal which 
has just two parameters: the interaction strength, /, controlled by the electron density; and 
the concentration of the minority spin state, p, set by the ratio of temperature to Zeeman 
splitting. It is similar to models for the band structure of uncorrelated binary alloys, and 
various approximate methods of solution are available. To facilitate a direct comparison 
with experiment, we choose, however, to solve Eq ^ by numerical diagonalization. 

Before presenting the results of our calculations, it is appropriate to review the ex- 
perimental data, the most striking aspect of which is the temperature dependence of the 
spectrum. At low density, two resonance peaks are observed at the positions of the separate 
spin transitions. Fig. 1(a). At a slightly higher density, the behaviour is much more com- 
plicated: both peaks shift continuously with temperature and eventually collapse together. 
Fig. 1(b). Above v ~ 1/6, there is a single peak (not shown), whose position depends weakly 
on temperature. 

These three types of behaviour correspond to the weak, intermediate and strong coupling 
regimes of our model. We have calculated the absorption spectra in each regime for a system 
of 225 particles, using periodic boundary conditions and Ewald summation for the dipole 
forces [|lT|, and averaging over 400 realizations of the spin distribution. Longer studies 



of up to 900 particles show that finite size and statistical errors are small under these 
conditions. We display our results in Fig 2. At weak coupling, the two spins are probed 
almost independently. Fig 2(a). At intermediate coupling, the resonance structure becomes 
strongly dependent on the spin polarisation, controlled by temperature. Fig. 2(b). In the 



strong coupling regime (not shown), we obtain a single resonance at the average cyclotron 
frequency, which shifts with temperature to follow the changing spin population. 

The successful reproduction in the calculated spectrum. Fig 2, of most of the observed 
features. Fig 1, constitutes our main result. We consider it persuasive evidence that our 
interpretation of the experiments incorporates the relevant physical ingredients. In making 
this comparison, the only scope for adjustment lies in the conversion between experimental 
electron densities and the interaction strength, J, of the model. The experimental densities 
of Fig 1 lead to values for the interaction strength of / = 0.1, 0.17, while we find the best fit 
with smaller values: / = 0.05, 0.1. We attribute the discrepancy to the finite thickness of the 
electron gas, which softens the Coulomb repulsion at short distances; comparable behaviour 
has been noted previously in calculations on the fractional quantum Hall effect [1^. The 



spin populations are determined by the Zeeman splitting, using the value [|T^, g = —0.4, for 
the electron (^-factor, and the theoretical spectrum is convolved with a Lorentzian to match 
the experimental resolution. 

It is initially surprising that the results of our calculations represent so well the observa- 
tions, since the theory assumes that guiding centre coordinates lie on a lattice, while most of 
the measurements are for temperatures well above the predicted melting point of the Wigner 
crystal. The explanation, we believe, is that the crossover from single-particle to collective 
cyclotron resonance is relatively insensitive to interparticle correlations, because the inter- 
actions are long-ranged. In principle, one could investigate resonance in the liquid phase by 
using molecular dynamics to follow guiding centre motion, simultaneously integrating Eqs ^ 
and §. The data presently available does not justify such an approach. Instead, we have 
examined the importance of interparticle correlations by diagonalising Eq ^ for fixed liquid 
configurations of guiding centres, generated using a Monte-Carlo method. This introduces 
a third parameter to the model, F = kT 47TeeQa/e^. In Fig. 3 we compare the results, at 
temperatures corresponding to 4K and 20K for an electron density of 3 x lO^'^cm"^, to those 
for the Wigner crystal. The lineshape changes by less than the experimental resolution, even 
at these high temperatures. 
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Further insight comes from studying a simphfied model (suggested by the long-range 
character of interactions) in which each electron is coupled equally to all others, with a 
coupling constant inversely proportional to the total number of particles. This model is 
obviously insensitive to the positions of guiding centres. It supports only two optically- 
active modes, the in-phase and out-of-phase oscillations of the two spin populations, with 
frequencies and oscillator strengths given by 



a;,„ = (a - 1 - ^/(l - «)' + 4ap)/2 (9) 



ooout = (« - 1 + V(l - «)' + 4ap)/2 (10) 

where p is the concentration of the minority spins. For comparison with our previous results, 
we choose the coupling constant, a, so that the sum of the interaction matrix over all 
neighbours is the same as that obtained with Coulomb interactions in the Wigner crystal, 
which requires that a = 11.034/. The frequencies and oscillator strengths are displayed in 
Fig 4: as the coupling is varied, there is a crossover from independent resonances of the 
two spin populations at weak coupling, to a single, in-phase, mode which will lie at the 
average cyclotron frequency, —p, at strong coupling. The crossover is most sharply defined 
for a single minority spin, p = 0, when it marks the transition from a bound state on 
the impurity spin (with relative oscillator strength inversely proportional to total particle 
number), to an extended state which carries all the oscillator strength. This transition 
would not be present in a crystal with short range interactions, since an arbitrarily weak 
attractive potential always has a bound state in two dimensions when the dispersion relation 
is quadratic at small wavevector. It does, however, occur in the two-dimensional Wigner 
crystal with Coulomb interactions, because the magnetoplasmon dispersion relation is linear 
at small wavevector. We find the critical interaction strength for the appearance of a bound 
state to be Ic = 0.104 in this case, while for infinite range interactions ac = 1, implying 
Jc = 0.091. 

Although spin-split cyclotron resonance is, as we have shown, very insensitive to spatial 
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correlations between electrons, it does provide a direct measure of spin polarisation. We 
demonstrate this by considering the first moment of the spectrum. To first order in 5uJc/oJc 



{w)= ujS{iu)diu/ S{iu)diu (12) 

~ Y^ {uJittia] + {l/h)diUa]) / J2(^i^l) (13) 

where di = (d/dxi+id/dyi) and U is the full interaction. For pairwise interactions, J2i 9iU = 
0, and in the absence of Landau level mixing, (aiOj) = 5ij, so we find 

1 ^ 
^^ i=l 

Thus the mean absorption frequency is simply the average cyclotron frequency of the parti- 
cles. Given the cyclotron frequencies for each spin orientation, it provides a direct measure 
of the relative populations. 

In conclusion, we have presented a model for spin-split cyclotron resonance at low filling 
fraction. Within this model we have been able to account fully for recent cyclotron resonance 
measurements on GaAs heterostructures. We find that the technique is rather insensitive to 
the liquid-solid phase transition, but provides a direct measure of spin polarisation. 

We are grateful to R. J. Nicholas for many helpful discussions. This work was supported 
in part by SERC grant Gr/Go 2727, and in part by DENI. 
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FIGURE CAPTIONS 

FIG. 1. Temperature variation of the experimental resonance spectrum, reproduced from 
[2]. The electron densities are (a) 2.6 x 10^°cm~^ and (b) 3.4 x 10^°cm"^. 

FIG. 2. Predicted resonance spectrum of the Wigner crystal, (a) / = 0.05, (b) / = 0.1 

FIG. 3. Comparison of the resonance spectrum in the Wigner crystal (F = 0) and in the 
liquid at two temperatures (F = 0.2, 1.0), for p = 0.1. 

FIG. 4. Crossover of the meanfield resonances, for p = 0, 0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5: 
(a) frequencies, (b) oscillator strengths. 
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